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PACS 05 . 70 . Ln - Non-equilibrium and irreversible thermodynamics 

PACS 05. 40. -a - Fluctuation phenomena, random processes, noise, and Brownian motion 
PACS 47 . 70 . Nd - Non-equilibrium processes, gas dynamics 

Abstract - An experimental study is presented, about transitions between Non-Equilibrium 
Steady States (NESS) in a dissipative medium. The core device is a small rotating blade that 
imposes cycles of increasing and decreasing forcings to a granular gas, shaken independently. 
The velocity of this blade is measured, subject to the transitions imposed by the periodic torque 
variation. 

The Hatano-Sasa (HS) equality, that generalises the second principle of thermodynamics to NESS, 
is verified with a high accuracy (a few 10 -3 ), at different variation rates. 

Besides, it is observed that the fluctuating velocity at fixed forcing follows a generalised Gumbel 
distribution. A rough evaluation of the mean free path in the granular gas suggests that it might 
be a correlated system, at least partially. 
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Introduction. — Recent decades have seen significant 
progress in nonequilibrium statistical mechanics, with the 
advent of the Fluctuation Theorems, the Jarzynski and 
Crooks relations [THS] . These relations were at the time 
theoretical advances, with the support of the numerics. 
Experimental contributions came later, mostly because of 
technical limits. Indeed, the scales at which thermal en- 
ergy dominates are small. Measurements of fluctuations at 
such scales has been prohibitively difficult until recently. 

Usually, an inequality involving the average entropy pro- 
duction is the expression of the second principle of thermo- 
dynamics. The improvement brought by the fluctuation 
theorems is that an instantaneous rate of entropy produc- 
tion is expressed by an equality. It is somehow a local 
formulation. 

The Jarzynski equality relates the Helmholtz free en- 
ergy difference AF between two states A and to the 
average of the exponentiated work needed to perform the 
transition: 



It can be equivalently written: 
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where Wdiss = W — AF is the work dissipated into heat 
during the transition. 

The Jarzynski relation is valid for any transformation, 
whatever the rate. For a reversible transition, Wdiss is 
obviously zero. The fluctuation theorems, as well as the 
Jarzynski and Crooks relations, refer to systems in equi- 
librium states, or submitted to transitions between equi- 
librium states, reversible or not. 

Another relation was derived latterly by Hatano and 
Sasa, in 2001. Generalising the Jarzynski equality, their 
prediction is drastically distinct as it addresses transitions 
between NESS of overdamped Langevin-type instead of 
equilibrium states. In that case, the forcing consists in 
nonconservative and potential forcings, together with a 
Gaussian white noise forcing, uncoupled to each other [4]. 
It writes similarly as the Jarzynski's equality (eq. |5J): 



(i) 
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The brackets denote the average over a large number of 
transition paths, and f3 = l/k^T, with &b the Boltzmann 
constant and T the temperature of the heat reservoir. 



with: 
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din [p ss (x; a) 
da 



(3) 



(4) 



p-1 



A. Mounier, A. Naert. 



The integral is evaluated over the transition time r be- 
tween two distinct NESS. The dot refers to time deriva- 
tive, and p ss (x;a) is the steady state probability density 
function (PDF) of the observable x at a specified value a 
of the control parameter. 

It is implicitly assumed that for any fixed value of a, 
the system relaxes to a single steady state characterised 
by p ss (x;a). Eq. [3] is expected whatever the transition 
rate. 

In a sense, the Jarzynski relation is a local extension of 
the 2 nd principle for equilibrium states, whereas Hatano- 
Sasa (HS) equation is its extension for NESS. Beyond the 
formal analogy between eq.[2]and[3j the HS relation comes 
from a distinct and more general phenomenological frame- 
work, called Steady State Thermodynamics [5j[6]. 

Regarding the prediction of eq. [3j only two experimental 
confirmations have been produced so far. First is that of 
Trepagnier et. at. [7] , that dragged periodically a colloidal 
particle in water, with an optic tweezer. This system ver- 
ifies all the requirements of the theorem, as the solvent is 
at equilibrium. It is a perfect case of Brownian motion, 
biased by an external conservative force. More innova- 
tive is the recent work of Gomez-Solano et. at. [8j. These 
authors performed a similar experiment, dragging a col- 
loidal particle with an optic tweezer in water. But instead 
of equilibrium states, it is prepared in NESS before cycling 
of the order parameter, according to a specified protocol. 

A step beyond, the present study is an experimental ev- 
idence that the HS equality also holds in a granular gas, 
i.e. for transitions between NESS in a dissipative medium. 

The HS relation (eq. [3|) refers to NESS Markovian pro- 
cesses, where fluctuations are not specifically of thermal 
origin. Therefore, the smallness of k^T must not be a 
limit... In other words, there is no need to study micro- 
scopic systems. The experiment presented here actually 
addresses a macroscopic system: it is extremely simple on 
its principle, and rather easy technically. 

In a dilute, continuously shaken granular g as, a blade 
is rotated around a vertical axis by a small DC motor 
at controlled torque. The angular velocity, resulting of 
this external torque and the numerous collisions with the 
beads, is the stationary fluctuating quantity under study. 
It is measured by the very same DC motor that forces 
the rotation. The torque, which is the control parameter, 
is ramped up and down periodically, causing transition 
between steady states of different mean velocities. His- 
tograms of the velocity are recorded for different values of 
the torque. 

These histograms appeared unexpectedly well fitted by 
a generalised Gumbel (GG) distribution. This is an in- 
termediate outcome of the present study, extremely useful 
for the calculations of eq. EJ Indeed, as the derivative of 
a GG distribution can be expressed exactly, the integral 
can be formulated. Therefore, it is easy to verify eq. [3j for 
different ramps of the control parameter. 

However, a tentative interpretation of this interesting 
observation is given in the last section. 



Experiment. — The set-up is sketched in fig. [TJ It 
is an improved version of the one used recently to study 
Fluctuation Theorem [9 . It makes use of a DC motor, 
converting current into torque, reversely used as a gener- 
ator to convert momentum into voltage. The same device 
is thus employed as actuator and sensor. A light plastic 
blade, embedded into a vibrated granular gas, is driven 
by a small and light DC motor. Forcing the rotation of 
the blade through the current (torque), one can measure 
simultaneously the rotation velocity through the voltage. 
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Fig. 1: The mechanical system is composed of a vibrating vessel 
containing the beads, excited by a shaker. The probing DC 
motor is fixed on the cover, here pulled out for clarity. 

The granular gas is composed of about 300 stainless 
steel beads of 3 mm diameter, vibrated in an aluminum 
vessel by a shaker. The vessel is 5 cm diameter and 6 cm 
deep, and its bottom is slightly cone-shaped to enhance 
horizontal momentum transfer. Thanks to a generator and 
a power amplifier, the shaker is supplied by a sine current 
at 40 Hz, providing a vertical acceleration of 41ms -2 . In 
that conditions, the granular gas is rather dilute. The 
blade is 2 cm x 2 cm, placed a few mm from the bottom. 
The nominal power of the DC motor is 0.75 W . The rotor 
is ironless, to minimise inertia, and precious metal brushes 
improve the electrical contact with the commutator. 

A current / injected into this motor results in a torque: 
T ex /, performing work against the granular gas. The 
same device can be used as a generator. In that case, 
the induced voltage e is proportional to the angular ve- 
locity: e oc 0. The proportionality factor accounts for the 
electro-mechanical characteristics of the motor. As it is 
the same in the motor or generator function, calibration 
is not needed. 

Note that the excitation of the vibrator that keeps the 
granular gas in a NESS by compensating the dissipation, 
is totally distinct of the torque applied by the motor to 
probe the gas. The former is a few Watts, as the later is 
a few mW to minimise perturbation as much as possible. 

The electric circuit is shown in fig. [2j The DC voltage 
supply is uo = 10 V (stabilised). The current / is driven 
by the voltage u supplied by a function generator. A time 
constant related to the inductance is irrelevant. The motor 
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Fig. 2: The electrical sketch of the motor's command. 

is depicted in fig. [2] as the assembly of a voltage source e 
(oc #), and the internal resistance r ~ 21.2 Q. The current 
is measured, thanks to a shunt resistor R = 56 Q. A 24 bits 
simultaneous data acquisition system records the signals 
1^0, u\ and ui at a sampling frequency of 1024 Hz. The 
instantaneous current I(t) and induced voltage e(t) are 
easily calculated from these voltage measurements: I(t) = 
(u (t) - uxit^/R and e(t) = u ± (t) - u 2 (t) -rl(i). 

The voltage generator is programmed to perform a cycle 
between low and high current regimes, i.e. torque cycles 
(fig. [3]). Various periods and transition rate have been 
performed, as discussed below. Each of the regime L or 
H, corresponds to a NESS. One precaution to be taken 
concerns the values of uo and u. They must allow currents 
/ such that the motor never completely stops rotating. 
Thus, no static friction is to be accounted for. Avoiding 
this difficulty is the reason why low values of torque are 
not explored in this work. 




Fig. 3: Sketch of the control parameter cycles, driving the sys- 
tem at 'High' torque and 'Low' torque NESS, or in constant 
rate transitions in-between. Transition time is r. 

A thermometer has been added on the vessel's cover 
to follow the temperature drift during the measurement. 
The temperature increases of about 5° during a typical 
transient time of 5 hours. This elevation of temperature 
perturbs the measurements, probably because of the vari- 
ation of air viscosity. Only the measurements performed 
after this transient of a few hours are considered. 



Principle. — When the torque is fixed, the blade ro- 
tates with a fluctuating angular velocity. The fluctuations 
are caused by the collisions with the granular gas. The 
equation of motion of the blade + rotor mobile writes: 



M6 + 1 {0)=T{t)+r 1 {t), 



(5) 



where is the angle, and dots stand for time derivatives. 
M and 7 are respectively the moment of inertia and the 
viscous friction term. This j{0) stands for an effective 
viscous damping due to the collisions, like in Brownian 
motion. It also includes a negligible (turbulent) drag on 
air. A small solid friction term is present, mostly in the 
commutator of the motor. It is constant, as the driving is 
such that the mobile never stops rotating. It can therefore 
be included as an offset in the torque T, and play no role. 
The deterministic torque T(t) is imposed from outside. 
The last term rj(t) is the random force accounting for the 
shocks of the beads. It represents the coupling with the 
NESS granular gas heat bath, i.e. the momentum transfer 
rate at each shock with the beads. All this description is 
written with Brownian motion theory in mind. Hints are 
given below that this description might not be correct. 

Eq. [5j that mimics the probing device, governs the ve- 
locity resulting from the balance between a deterministic 
forcing, the coupling with a stady state reservoir, and 
friction. At first glance, it takes the form of a Langevin 
equation, if the noise 77 can be considered short-time 
correlated. However, a first difficulty comes from the 
dependences of the forcings with one another. Indeed, 
the random force 77(f) is affected when T(t) is changed, as 
shown below. The whole balance between deterministic 
and random forcings is varied. The angular velocity 
follows in a non trivial manner. All things considered, 
the description of this system with eq. [5] as a Langevin 
equation is not as simple as it first appears. 

The purpose of this work is to study transitions 
between two NESS, characterised by the fluctuating 
angular velocity 6(t), while T(t) is ramped at fixed rate 
between two specified values corresponding to 'states' 
L and H. The HS equality is verified with a very good 
accuracy. 

For convenience, another set of variables than T} 
is used. The observable e (in Volts) is centered and 
normalised, such as x(t) = (e(t) — e)/cr, with the mean 
e, and the variance a 2 = (e — e) 2 . (The bar denotes 
time- average within a single steady state.) The control 
parameter is from now on the current I(t) (in Amperes). 

As already mentioned, the calibration factor is left 
aside, not necessary to test eq. [3) This equation is 
rewritten with the new working electric variables: 



exp 



dtl 



•din \p aa (x]I)] 



dl 



1. 



(6) 



As discussed above, the quantities {x,I} are directly mea- 
sured. The whole analysis procedure is performed on this 
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new set of variables. The derivative of In [p S s( x ; I)] is to be 
taken from the histograms calculated over large samples 
of x, at fixed values of /. 

The integration over the transition time is easy, as well 
as the average over a large number of transitions. The 
derivative of histograms with respect to the control pa- 
rameter / is actually more difficult. However, a specific 
character of these histograms is to be well fitted by a GG 
distribution (see next section). This observation is of great 
help for the analysis procedure. 

The generalised Gumbel distribution. — The his- 
tograms of x at fixed / have an asymmetric but univer- 
sal shape, whatever the value of the control parameter /. 
They are very well fitted by a GG distribution (fig. [4]). 




2 4 

(e-e)/a 

Fig. 4: A histogram of the centered and normalised induced 
voltage x(t) — (e(t) — e)/a at fixed current / is plotted in semi- 
log axis (dots). The fitting is performed with a GG distribution 
(line). The best fit is obtained for a — 2.5. 

Assuming the variable x is distributed according to a 
GG law, it is characterised by a single shape parameter a, 
that accounts for the asymmetry: a ~ l/< x 3 > 2 . (The 
PDF tends to a Gaussian distribution if a —> oo.) It writes: 

p ss (x) = K a exp [a [-b a (x + s a ) - exp (-b a (x + s a ))}} , 

(7) 

The mean and the variance, as well as the normalisation 
factor, can all be expressed as functions of a: 



ba 



d 2 lnr(a) 



da 2 



l/, n dlnr(a) 
Sa = T ln ( a ) a 
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(8a) 
(8b) 
(8c) 



thanks to the gamma- function: T(a) = f™t a - 1 e~ t dt. 

As the PDF's shape as well as the mean and stan- 
dard deviation only depends on a, it can be rewritten as 
p ss (x;a) after a change of variable. The integral term Y 



of eq. [6] is therefore rewritten 
Y -- 



dt I ( — \ d ln [pss(x] a)] 
\dl J da 



(9) 



The dependance in / of the parameter a, obtained from 
the fitting of the histograms, allows to calculate (^f). 

Assuming p ss is a GG distribution, the differentiation 
needed in eq. [9] can be performed exactly. 

Results. — To center and normalise the variable e, the 
mean e and standard deviation a are directly measured 
from the voltage time series corresponding to stationary 
states for over ten fixed values of /. They are plotted 
against the current / in fig. [5] A linear fitting is performed, 
valid at least in the range of interest. 
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Fig. 5: The mean value e (up) and the standard deviation a 
(down) at fixed current / is plotted against /. A linear fitting 
is performed over the range available. 

The standard-deviation a is a growing function of /, not 
expected to cancel for / = 0. Indeed, at zero-torque, fluc- 
tuations of velocity remain, because of the random forcing 
r\. By effecting certain calibration, it could be linked in a 
non local manner to a granular temperature [TO] . 

It is to be noticed that an extrapolation of e following 
the linear fitting does not go to for I = 0. This is ob- 
viously abnormal, as the blade should not rotate without 
torque (for symmetry reason, as < rj >= 0). There might 
be a nonlinearity j(0) in the 'viscous drag' of eq. [5] This 
point is discussed below. 

Now, the fitting of the histograms for different values 
of I is performed, and the parameter a is extracted. It is 
plotted against / in fig. [6j This parameter a increases for 
lower /, meaning that the distribution symmetrises when 
the external excitation decreases. Joubaud et. al. recently 
observed in a granular gas, that velocity fluctuations with- 
out external forcing look Gaussian [11] . Their experiment 
is designed for small static friction. It is not quite ac- 
curate to do such measurements here in this low-torque 
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regime, because of the friction in the commutator of the 
DC motor. 



Table 1: The HS eq.[3]is confirmed with a very good accuracy: 
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Fig. 6: The asymmetry coefficient a is obtained for each I 
by fitting the histogram with a GG distribution, a is plot- 
ted against the current 7, together with the inverse squared 
skewness. 

The linear fittings give the simplest dependance of those 
three quantities with 7: 



e = 
a 



9.3 7- 
-2.61 



0.12, 
-3.2 1(T 
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It is not surprising to notice that the statistical noise is 
larger for increasing order moments: e, a, and a. 

A tedious derivation leads to the following exact expres- 
sion for the log-derivative of the GG distribution: 

dln[p ss (x;I)] _ 



da 



1 

2¥ 



_ 2x ^ 2 + (a - e^- x ^)(2^ /2 - xV^ 7 *") 



(11) 



where \I/(a) = dlr ^ a ^ is the so-called digamma function, 
and ^" its successive derivatives with respect to a. 

In the limit a — >• oo, this expressions reduces to a 
quadratic form: din[ps S jx,i)] _ _L(£ 2 — 1), consistent with 
a Gaussian p ss (x;I), as < 0. This limit refers to van- 
ishing 7, as discussed above; the linear relation between 7 
and a is only a working approximation. 

The log-derivative is computed for all the measured time 
series x(t), then multiplied by I and ^j. Realisations of 
Y are obtained by integration for each transient. Thence, 
the average of the exponential over dozens of transitions is 
carried out, separately for leading and trailing edges (in- 
creasing and decreasing torques), for slow and fast transi- 
tion rate. Results are shown in the following table: 





leading edge 


trailing edge 


r = 10 s 


0.9890 


1.0069 


r = 30 s 


1.0012 


0.9985 



Y is equivalent for one path between two NESS to the 
dissipated work between two equilibrium states calculated 
through Jarzynski equation. It would certainly be inter- 
esting to compute its histogram. However, the sample 
would have to be much larger than that presently avail- 
able. 

Discussion. — This article presents an experimental 
study of a granular gas, regarded as an ersatz of a heat 
reservoir. The granular gas is considered as a thermostat, 
however dissipative. A simple device coupled to this reser- 
voir exchanges energy with it. 

This experiment takes advantage of the fact that small- 
ness of the systems is not required. The granular gas is 
probed with a blade rotating about its vertical axis, which 
velocity is measured at controlled torque. The torque is 
cycled in such a way that angular velocity undergoes tran- 
sitions between stationary states. 

The Clausius inequality gives a lower bound to the work 
dissipated in transitions between equilibrium states. The 
HS equality generalises it to transitions between NESS. 
This article describes the first experimental observation of 
the HS prediction in a dissipative system. The agreement 
is impressive even if the dependence of the distribution is 
not known over the full range. It appears to hold indiffer- 
ently whether the forcing is undergoing an increasing or 
decreasing transient, whether steep or gentle. 

Strictly speaking, the HS relation is expected to be valid 
for stationary states. However, one could expect no de- 
parture as long as the transition time r is larger than a 
microscopic time of the reservoir's fluctuations, where re- 
arrangements can occur during the evolution of the order 
parameter. As it is the mean time between two shocks, 
such rapid transition is probably limited in the present 
experiment by the inertia of the blade. Therefore, the 
range of applicability of the HS relation is larger than ex- 
pected, from this point of view. A generalisation of HS 
theorem to non stationary processes is discussed in [12] . 

It is assumed from the beginning that the gas is dilute. 
To make this statement quantitative, the mean free path is 
evaluated, thanks to crude dimensional arguments. First, 
the density. Because of vertical stratification, density is 
larger in the lower part of the cell. If all the beads are 
assumed uniformly distributed in the lower h = 1 cm of 
the cell, the density is n = N/(7rR 2 h), with R the radius 
of the vessel and N = 300 the number of beads. It gives 
n ~ 16 cm -3 , which means that the mean distance be- 
tween beads is about 4 mm. Now, following the kinetic 
theory of gases, the mean free path is: A = l/(n7rr 2 ), 
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where r = 1.5 mm is the radius of a bead. It gives 
A ~ 9 mm. This evaluation is a rough order of magni- 
tude, it should be improved. However, it is close to any 
length in this experiment! 

Besides, the correlation time of x(t) is of the order of 

24 ms. This corresponds roughly to \/{L6). Correla- 
tion in space and time are consistent. 

This result, associated with the asymmetry of the PDF 
and the nonlinearity of the 'drag' 7(0), shows that a de- 
scription of this system in terms of a simple Brownian 
motion is too simplistic. The Knudsen number must be 
considered, defined as: Kn = A/L, L being a characteris- 
tic length of the system, like the radius of the blade. 

The central requirement for a process to verify the HS 
equality is to follow an overdamped Langevin equation. 
The gas being rarefied means that collisions come one by 
one on the blade. Therefore, the noise n(t) is not likely 
to be a Gaussian white noise, and assuming a (nonlinear) 
friction j(x) is not representative of the physical reality. 
There is no reason to neglect the inertial term in eq. [5j For 
all these reasons the equation of motion does probably not 
verify required conditions, i.e. with such a high Kn, the 
system does not behave like simple Brownian motion. 
In such case, the experimental verification of the HS equal- 
ity enlarges its range of validity. A better determination 
of A, or a direct test of the Markovian character of the pro- 
cess is badly needed. (Numerically?) An opening would 
be trying this relation with experimental processes clearly 
non- Markovian, or rapidly varying, change parameters like 
stratification, density, sizes, excitation, to identify which 
conditions causes failure of HS prediction. 

Besides, it is shown that the fluctuations of velocity at 
fixed forcing, and therefore the power injected by the blade 
into the granular gas, are asymmetric and resemble a GG 
distribution. Such kind of distribution have been found 
describing fluctuations of power injected in dissipative or 
correlated systems (see experiments on turbulent flows in 
[T3l[l4], numerical simulations on granular gases in [15]), 
or other global quantities such as the fluctuations of mag- 
netisation in critical ferromagnetic systems with finite size 
effects (see XY or Ising model computations in [T6| [T7] ) . 
The usual theoretical explanation for such statistics is that 
global quantities's fluctuations are affected by correlations 
(see [18j[T9] and references therein). The very basic idea 
is that the shape parameter a is linked to the number of 
degrees of freedom of the system: asymmetry comes from 
the finiteness of this number. In the present situation, the 
tentative explanation rely on a dependance between the 
relatively large value of Kn and correlation. 

Another cause of asymmetry in statistics can be clus- 
tering, due to dissipation. However, it would simply en- 
hance vertical density stratification, without causing clus- 
ters that the blade would hit during its rotation. 

A blade of half width (1cm x 2 cm) has been tried in 
the same configuration. As a result, the fluctuations are 
much more asymmetric (a is much smaller, for instance 



5.52 — )• 2.45). This observation corroborates qualitatively 
the previous argument, as Kn is doubled. 

At this point, it is important to clarify the Markovian 
character of the process involved, and the correlations in 
the system. Answers to these questions could mean a 
widening of the conditions of this theorem, and explain 
the asymmetric statistics altogether. It could also be in- 
teresting to relate the parameter a to Kn. 

* * * 
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